Propagation effect on polarization of pulsar radio emission 



R. T. Gangadhara 1 ' 3 , H. Lesch and V. Krishan 1 

1 Indian Institute of Astrophysics, Bangalore-560034, India 

2 Institut fur Astronomie und Astrophysik der Universitdt Miinchen, Scheinerstrasse 1, 81679 Miinchen, Germany 
3 National Centre for Radio Astrophysics, TIFR, Pune University Campus, Pune-J^llOOl , India 



Received: 

ABSTRACT 

We consider the role of a propagation effect such as the stimulated Raman scattering on 
the polarization of radio pulses. When an intense electromagnetic wave with frequency 
close to the plasma frequency interacts with the plasma in the pulsar magnetosphere, 
the incident wave undergoes stimulated Raman scattering. Using typical plasma and 
magnetic field parameters, we compute the growth rate and estimate the polarization 
properties of the scattered mode. At some conditions, we find that the polarization 
properties of the scattered mode can become significantly different from those of the 
incident wave. The frequencies, at which strong Raman scattering occurs in the outer 
parts of magnetosphere, fall into the observed radio band. In some pulsars, for example, 
PSR B0628-28 and 1914+13, a large and symmetric type of circular polarization has 
been observed. We propose that such an unusual circular polarization is produced by 
the propagation effects. 
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1 INTRODUCTION 

The investigation of coherent radio emission from the pulsar magnetosphere has attracted a great deal of attention (Cordes 
1979; Michel 1982; Asseo et al. 1990). A powerful collective emission occurs when relativistic electron beams with density ~ 
1 per cent of the pair plasma density, scatter off coherently from concentrations of plasma waves (cavitons) (Benford 1992). 
The role of collective plasma processes in the absorption and spectral modification of the radio waves is well known (e.g., Beal 
1990; Krishan & Wiita 1990; Benford 1992; Gangadhara & Krishan 1992, 1993, 1995; Gangadhara, Krishan & Shukla, 1993). 
von Hoensbroech, Lesch and Kunzl (1998) have demonstrated that degree of linear polarization decreases with increasing 
frequency while the degree of circular polarization shows the opposite trend. 

In this paper, we estimate the role of stimulated Raman scattering on the polarization of electromagnetic waves prop- 
agating in the pulsar magnetosphere. We assume that the physical conditions in the pulsar magnetosphere are those of the 
classical standard model (Ruderman & Sutherland 1975) which describes the generation of ultrarelatvistic beams of elec- 
trons/positrons and the creation of the pair plasma. The beams and the pair plasma are in relativistic motion along the 
bundle of open magnetic field lines that delimit the active region of magnetosphere. The stimulated Raman scattering is 
considered as a parametric decay of the initial transverse electromagnetic (pump) wave into another electromagnetic wave 
and a longitudinal plasma wave. The physics of stimulated Raman scattering in a plasma has been explained in many papers 
and books (e.g., Drake et al. 1974; Liu & Kaw 1976; Hasegawa 1978; Kruer 1988). In two ways, stimulated Raman scattering 
may be important in the pulsar environment: first, it may act as an effective damping mechanism for the electromagnetic 
waves generated by some emission mechanism at the lower altitudes in the pulsar magnetosphere. At those altitudes, the res- 
onant conditions for stimulated Raman scattering, i.e., the frequency and wave number matching might not be satisfied. This 
results in a short time variability which is generally observed in pulsar radio emission. Secondly, it may provide an effective 
saturation mechanism for the growth of the electromagnetic waves provided that the conditions for the wave excitation by 
some mechanism are satisfied in the region where an effective stimulated Raman scattering can take place. 

The first case can be simplified by treating the intensity of the pump is constant in time. Then the nonlinear equations, 
which describe the wave coupling, become linear in amplitudes of the decaying waves, and the exponentially growing solutions 
will imply an effective energy transfer from the pump wave. This approximation breaks down when the amplitudes of the 
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decay waves become comparable to the pump wave or when the amplitudes of the decay waves enter nonlinear stage, and 
start loosing energy due to some nonlinear processes such as wave-particle trapping and acceleration. 

The second case is more complicated, where stimulated Raman scattering acts as a nonlinear saturation mechanism and 
the amplitudes of all waves may be of the same order. This case can be considerably simplified when damping of the plasma 
wave is very strong or if it leaves the region of the resonant interaction fast enough. 

We neglect the nonlinear stages of stimulated Raman scattering and the development of Langmuir turbulence, which 
leads to wave-particle trapping or quasilinear diffusion. If the pump is monochromatic, the growth rate of stimulated Raman 
scattering can become very high, as in conventional laboratory laser-plasma interaction. However, in the case of pulsars the 
pump can be broadband, and in the limit where the bandwidth Aui of the pump wave is much larger than the growth rate 
of stimulated Raman scattering, we can use a random phase approximation for the statistical description of the interacting 
waves. 

Tsytovich and Shvartsburg (1966) have given a general expression for the third order nonlinear currents excited in a 
magnetized plasma. Since the corresponding expressions are very complicated, the general case of Raman scattering becomes 
very difficult to consider. However, one can make some useful simplifications, when considering stimulated Raman scattering in 
the pulsar magnetosphere. First, in superstrong magnetic field, we can expand the currents in 1/lub, where lob is the cyclotron 
frequency. Second, if the pair plasma has the same distributions for electrons and positrons then some of the nonlinear currents 
cancel out, as they are proportional to the third power of the electric charge (this cancellation is exact in the unmagnetized 
electron-positron plasma). Third, all the three interacting waves propagate along magnetic field. This is an important but less 
justified approximation. It allows us to simplify the problem considerably, and to obtain a dispersion relation for stimulated 
Raman scattering. 

The polarization of the pulsar signals appear to bear critically on the pulsar radio emission process and the emission beam 
geometry. One or more reversals of the sense of circular polarization has been observed in the intergrated profiles of several 
pulsars. However, in individual pulses circular polarization changes sense many times across the pulse window. Further, it is 
important to determine whether the depolarization is a geometric effect or results from radiation-plasma interactions. There 
have only been very preliminary attempts to explain depolarization and microvariability using plasma mechanisms (Benford 
1992). 

Our purpose in this paper is to derive the properties of natural plasma modes and to explore some possible implications 
concerning the interpretation of the observed polarization, notably large and symmetric circular polarization in some pulsars 
(e.g., PSR B0628-28, and 1914+13). In §2, we derive the dispersion relation for stimulated Raman scattering of an electro- 
magnetic wave, and give an analytical expression for the growth rate of the instability. In §3, we define the Stokes parameters 
and compute the polarization states of the scattered electromagnetic wave. The discussion of our findings is given in §4. 



2 POLARIZATION CHANGES DUE TO STIMULATED RAMAN SCATTERING 

We begin with a model consisting of a pulsar with nonthermal component of radiation interacting with plasma in the emission 
region at a distance r — WORns = 10 8 cm, (neutron star radius Rns ~ 10 km) from the neutron star surface, where the 
magnetic field is about 10 6 G. Plasma particles may be all in their lowest Landau level with no Larmor gyration, however, 
plasma can have one-dimensional distributions of momenta along the magnetic field (Blandford 1975; Lominadze, Machabelli 
& Usov 1983). This is because the synchrotron loss time for the decay of the perpendicular component of momentum is 
always short compared to the transit time at the stellar surface for any velocity. The simplest model for the plasma is density 
declining in proportion to r -3 with no gradients in the distribution functions across B. 

The nonlinear interaction of radiation with plasma causes the modulation instability leading to enhancement of nonres- 
onant density perturbations and the radiation amplification by free-electron-maser, which produces intense electromagnetic 
waves (Lesch, Gil & Shukla 1994; Gangadhara, Krishan & Shukla 1993; Gangadhara & Krishan 1992). Since the frequency 
of these electromagnetic waves is close to the plasma frequency, they resonantly couple with the subpulse-associated plasma 
column in the pulsar magnetosphere and drive the stimulated Raman scattering. 

Consider a large amplitude electromagnetic wave (ki,u>i) with an electric field: 

Ei — E X i cos(fo • r — unt)e x + £ y i cos(fo • f — u>it + 8i)e y (1) 

interacts with the plasma in the pulsar magnetosphere. 

We follow Ruderman & Sutherland's (1975) approach to estimate the density and plasma frequency of the plasma moving 
within the bundle of field lines. For the typical parameters: Lorentz factor of primary particles 7 P ~ 10 , and for pair plasma 
particles 7± ~ 10 3 , magnetic field B ~ 10 6 r^~ 3 G, and pulsar period p — 1 s, we obtain the particle number density 

7p B 3.5 x 10 7 _ 3 

n = = s cm (2) 

2-y± ecp ?g 

and the plasma frequency 
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Figure 1. Stimulated Raman scattering of transverse electromagnetic wave (E x i, E y i) through an angle <j> a . The scattered wave electric 
field is (E xs , Eys)- The wave numbers ki, k a and k are due to the incident, scattered and Langmuir waves, respectively. 



w p = 27± [ ] = 3^ — rad s . (3) 



7±m 



The plasma in the emission region of pulsar magnetosphere may be birefringent (Melrose & Stoneham 1977; Melrose 1979; 
Barnard & Arons 1986; von Hoensbroech, Lesch and Kunzl 1998). In these models, two modes of wave propagation is generally 
allowed in magnetoactive plasma: one mode is polarized in the ki-B plane and other mode in the direction perpendicular to 
it. Following these models, we resolve Ei into two modes E xi — e xi cos(fc 4 • f — u>it)e x and E y i = e y i cos(fc 4 • f — u>it + Si)e y 
such that they are polarized in the directions parallel and perpendicular to the ki-B plane, respectively. It is well known from 
laser-plasma interactions that large amplitude electromagnetic waves resonantly interact with the plasma when the radiation 
frequency becomes close to the plasma frequency. Since the two modes have different indexes of refraction (McKinnon 1997), 
i.e., the response of plasma is not same for the two radiation modes, it is reasonable to assume that E xi couples with the 
density perturbation Am = 5ni cos(fc • r — cot), and E y i couples with A712 = 5n2 cos(fc • r — ut + S) in the plasma medium. 

Since the ponderomotive force is proportional to VE^ and VEf 2 , the coupling between the radiation and the density 
perturbations is nonlinear. Hence the density perturbations grow up which lead to currents and mixed electromagnetic- 
electrostatic side-band modes at (ki ± k, u>i ±u>). In turn these side-band modes couple with the incident wave field, producing 
a much stronger ponderomotive force, which amplifies the original density perturbation. Hence a positive feedback system 
sets in, which leads to an instability. 

The electric field E s of the electromagnetic wave scattered through an angle <j> s with respect to ki can be written as 

E a = e xs cos(fc s ■ r — cu s t)e x + e ys cos(fc 3 • r — u) 3 t + S s )e' y . (4) 

The propagation directions of incident wave {ki,u)i) and scattered wave (k s ,u> s ) are illustrated in Fig. 1, such that ki || e z , 
k s || e' z and e' y \\ e y . The primed coordinate system is rotated through an angle 4>s about the y-axis. Then the scattered wave 
in the unprimed coordinate system is given by 

E s = e xs cos(0 s ) cos(fc s ■ r — LJ s t)e x + e ys cos(fc s • r — uj 3 t + 8 s )e y — e zs sin(0 s ) cos(fc s • f— uj 3 t)e z , (5) 
where e zs = £ ls . 

The quiver velocity u± of positrons and electrons due to the radiation fields Ei and E s is given by 

^ = ± ^ + ^ ^ 
where e and mo are the charge and rest mass of the plasma particle. 

The wave equation for the scattered electromagnetic wave is given by 

where c is the velocity of light and J is the current density. 
The components of equation (7) are 

D s e XB cos(0 s ) cos(fc s • r — uj s t) = —^^——e X i8ni cos(fc s • r — w s t) (8) 

mo u)i 

and 
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D s e ys cos(fc s • r — uj s t + S s ) = -£ y i8ri2 cos(fc s • r — io s t + Si — 8) (9) 

mo Ldi 

and 

D s e zs sin(0 s ) cos(fc s • r — u a t) = 0, (10) 
where D s = k 2 s c 2 — uj 2 + ui 2 and 

lo s — LUi — ui, k s — ki — k. (11) 

In quantum language these two relations may be interpreted as the conservation ol energy and the momentum along the 
magnetic field, respectively. When these conditions are satisfied, stimulated Raman scattering is exicted resonantly and the 
expression D s w becomes the dispersion relation of the scattered electromagnetic mode. 

If we cancel the instantaneous space and time dependent cosine functions on both sides of equation (8), we get 

D 3 e xs cos<b s = -e x i8n\. (12) 

Similar to the phase matching conditions (equation 11), equation (9) gives the condition between the initial phases: 

S s = S t - 5, (13) 

and hence we obtain 
2?re 2 tj s 

D s e ys = e y idn 2 . (14) 

mo u>i 

Dividing equation (14) by equation (12), we have 

q s = ai7)cos(0 s ), (15) 

where a< = e y i/e X i, a s = e ys /£xs and ij = 5ri2/5ni. The value of r\ is determined by Ei in such way that Ani couples with 
E x i and Ari2 couples with E y i. 

We consider the Vlasov equation to find the low frequency plasma response: 

K-+^.V/ + — (eV0-V^).f£ = O, (16) 
at mo ov 

where <f>(r, t) is the scalar potential associated with the electrostatic waves, f(f, v, t) is the particle distribution function and 
ip(f,t) is the ponderomotive potential. 

Using f(r,v,t) = fo(v) + A/i(r,w,t) + Af2(r,v,t), we can linearize equation (16), and obtain: 

°mi + ?im + , v(i/l , + „,v,a/ 5 , + ±w - v*.f . o. m 

where A/i = 5/i cos(fc-r — ut) and A/2 = 5/2 cos(fc-r — wt+5). The ponderomotive force of the radiation field depends quadrat- 
ically on the amplitude of the electric field. Physically, it is a radiation pressure which amplifies the density perturbations by 
exciting the slow longitudinal fields and motions. The ponderomotive potential is given by: 



* = 2^\l Re 



= 7 [cos((f>s) cos(fc • r — uit) + aids cos(fc • r — cut + Si — S s )]eiS 3 , (18) 

2mouJi 

where the bracket ^ represents the average over the fast time scale (wj 3> oS). 

In the presence of a strong magnetic field, the pair plasma becomes polarized and one dimensional. If there is some 
relativistic drift between the electrons and the positrons then the ponderomotive force becomes effective, which will induce 
the nonlinear density perturbations (Asseo 1993): 

6H± 327rjk B T pW ?7±' 1 ' 

where eu is the envelope of the parallel electric field of Langmuir wave which is slowly varying with space and time, and T p is 
the plasma temperature. Using these density perturbations and Poisson equation, we can self-consistently derive the potential 

= -^(Ani + An 2 ). (20) 
Now, substituting the expressions for A/i, A/2, Am, Ari2, 4> and V 1 m to equation (17) and using the condition (13), we obtain 
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5f 2 + 5fm + 



4ire e 

— fe2 (5n 2 + 5ni/z) + r,^ 2 , ,2 {cos(<ft 3 )^ + aia s }eie s 



2m 2 u; 2 



k 9 f° 

dv 

(uj — k ■ v) 



0, 



(21) 



where /_t = sin(fc • r — u>t)/sm(k ■ f— cut + S). For 5 = and tv, fi — ±1, while at the other valuse of S we have to find the 
average < /j >= fi a over the time scale T = 2-k/u. Therefore, we have 



Sf'2 + SflfJ, a = - 



Aire 2 

7710 fc 2 



5U2 + Sn!fl a + 



Sivmocof 



k a ?5 

dv 



(lo — k ■ v) 



where A = (fj, a cos(f> s + aia s )e s . The sum of density perturbations (Sri2 + 5nifj, a ) can be estimated as 

oc 

Sn 2 + SniHa = J no(Sf 2 + 5fiiJ, a )dv = 



8n 2 + Sn-i^a + - £lk — o A 
87rmot<j 



X, 



where 



X ~ fc 2 



k 9fo 

— dv, 



(lo — k ■ v) 



is the susceptibility function (Liu & Kaw 1976; Fried & Conte 1961). Using r\ = 5ri2/5ni, we can write 

Sik i 



1 + — ) (v + Va)Sni 



8irniouj 2 

Multiplying equation (12) by Eifx a and equation (14) by ctiSi, and adding we get 



(a s Qi + Ha cos 4>s)e s 



2-Ke 2 lo s 2 sSm 

ei(Ha + ai7)) — , 

m UJi Us 



where Ei = e x % and e s = e xs - 

Now, using equations (25) and (26), we obtain the dispersion relation for stimulated Raman scattering: 

t > 2 fc 2 (fig + a\rj) u) s 1 
4 (1 + Q 2 ) UiDs' 



1 + - ](r? + At a ) = 

A, 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



where vo = eEi y/l + a 2 / rriQUJi is the quiver velocity of plasma particles due to the electric field of incident electromagnetic 
wave. 

If L = L30 x 10 30 erg s _1 is the luminosity of pump radiation with frequency uii = 27rz/» = uio x 10 10 rad s" 1 at a distance 
r = rs x 10 8 cm from the source then 



wo = 



mouii \ r c 



2L 



1/2 



= 4.3 x 10 



7-1/2 

9 i^- cms" 1 . 
triors 



(28) 



When the phase matching conditions are met, the instability becomes more efficient, and D 3 = 2a>;u; — w — 2c 2 fci.fc + c 2 fe 2 « 
represent the dispersion relation of the Stokes mode. 

The thermal speed w t of the plasma can be expressed in terms of its energy spread in the laboratory frame. Let v± = 
c(l— 1 /7± )*' 2 be the velocity of electron-positron plasma, then the velocity spread 5v± is given by (Hasegawa 1978; Gangadhara 
& Krishan 1992) 



for 7 > 1. 



Now, using the Lorentz transformation of velocity, we can show that 
8v' z _ <5«^ « t 



5u± = <5w z = 



7|(l + « ± ^/c 2 ) 7 2 



because v' z — 0. Hence the thermal speed, in the plasma frame, is given by 
#7± 



Vt = c- 



7± 



(29) 



(30) 



(31) 



For 7± = 10 3 , we get » t = 3 x 10 7 57± cm s _1 . 

If we separate equation (27) into real and imaginary parts, we get two coupled simultaneous equations. By solving them 
numerically, we find the growth rate T of stimulated Raman scattering. Figure 2 shows the growth rate as function of rg 
and uii/uip. Since the plasma density decreases with r as 1/r 3 , growth rate decreases as r$ increases. Also, if the radiation 
frequency becomes higher than the plasma frequency, then radiation and plasma do not couple resonantly, which leads to the 
decrease of growth rate at higher values of cd; /ui p . 
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Figure 2. The growth rate T of stimulated Raman scattering is plotted with respect to r$ and u>i/u) p . 

To study the variation of F in the ai-5'y± plane, we have made a contour plot (Fig. 3), the labels on the contours indicate 
the values of log(r). Debye length increases with the increase of ^7±, therefore, plasma wave experiences a large Landau 
damping, which leads to the drop in growth rate. The parameter a; is the ratio of amplitudes of electromagnetic waves, 
which are polarized in the directions parallel and perpendicular to the fc; — B plane. When a; is small, the wave polarized in 
the direction parallel to ki — B plane becomes strong and will have a component along B. Therefore, the coupling between 
radiation and plasma will be strong, which leads to the higher growth rate at the smaller values of a;. 
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Figure 3. Contour plot of T in the S^y±-ai plane. The labels on the contours indicate the log(T) values. 

For u> c 2 ki ■ k/oJi, Stokes mode becomes resonant, and anti-Stokes mode becomes non-resonant. Then the equation (27) 
can be written as 
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(a; - un + iT,)(a; - a,, + iT.)fo + Ma) = -10 18 7 ff ^ ± If rad 2 s~ 2 , (32) 

r 8 ' i^io U + ^J 



where ujf = u p + (3/2)fc i>f , and 



r ' = — 7T-; n 3 ex P 



2 (fcA D ; 



3 



2(fcA D ) 2 2 



+ (33) 



is the damping rate of the plasma wave, and the Debye length \q = v t /^/2tu p . For k\o ~ 1, we find T; « 2.5 x f0 9 /r| 2 s _1 . 
The pair plasma collision frequency v c « 2.5 x 1CP 3 In A/r|<57± s _1 and Coulomb logarithm In A w fO. The collisional damping 
rate of the scattered electromagnetic wave is given by T s = ujpV c /2iu1 w 0.06/f? ^8^7± s _1 . 

For fcAo > 1, plasma waves are highly damped and hence plasma losses its collective behaviour. Therefore, stimulated 
Raman scattering changes into the induced Compton scattering of electromagnetic waves by the plasma particles. In the 
conventional treatment of induced Compton scattering in pulsars (e.g., Blandford & Scharlemann 1976; Sincell & Krolik 
1992) the collective effects of the plasma are ignored. The collective treatment of the wave scattering by plasma particles is 
justified if the condition k\D « 1 or v p h >• vt is met, where v p h = w p /k is the phase velocity of the plasma wave. This 
condition implies that the wave number of oscillation of the electrons in the beat wave of the incident and scattered waves be 
much less than the inverse of the Debye length. It is not satisfied, when the beat wave is strongly Landu damped (k\o > 1) 
or the beat wave will not feel the presence of a medium (k\o ^ 1), so that the scattering process will be described as induced 
Compton scattering. 

When k\o <C 1, by setting u> = loi + iT, we can solve equation (32) for the growth rate: 



+r -' ± W (r ' - T ->' +43 * 10 " (i i%£t.y <34) 



Stimulated Raman scattering is a threshold process: if the intensity of the pump exceeds the threshold, then only it would 
start converting its energy into the decay waves. The threshold condition for the excitation of stimulated Raman scattering 
is given by 

(^) =9 .3xlO- 1 Vor 8 3/2 r,r a (1 + a?)(,? 2 + r ) - (35) 

The typical threshold intensities for stimulated Raman scattering are of the order of the observed intensities, implying that 
pulsar magnetosphere may be optically thick to Raman scattering of electromagnetic waves. 
The growth rate just above the threshold is given by 

r = 4.3xl0 8 L 2 30 - iEl+^Hl s ~\ (36) 
rjjvio (l + a 2 )(7? + Ma) 

which is proportional to L30. The maximum growth rate attainable for uj p > Y > Fi, on the other hand, is 



r-10 9 / L 3o (/jq + Q^) -1 r37) 



3 STOKES PARAMETERS 

When the phase matching conditions are satisfied, the growth rate T becomes large, and the scattered mode is amplified and 
become a normal mode of the plasma. Under some conditions, the scattered mode leaves the plasma with polarization which 
may be different from the polarization of the incident wave. The polarization states of the incident and scattered waves can 
be described more accurately using the Stokes parameters (Rybicki & Lightman 1979): 

Ij = £xj£ x j + e yj£yj> (38) 
Qj — £xj£ x j ~ £ yj£yj> (39) 
Uj = 2e x j£* yj cos{5j) (40) 
and 

Vj = 2e xj el j sm(S j ), (41) 
where j = i for the incident wave, and s for the scattered wave. The linear polarization is given by 



U'f + Ql (42) 
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and the polarization position angle is given by 



X j = iarctan^/Q,). (43) 

The transfer of energy between the modes will be efficient only when the energy of the pump wave is strong enough to 
overcome the damping losses or escape of the generated waves. Using the Manley-Rowe relation (Weiland & Wilhelmsson 
1977) 

- = -, (44) 
we find the relation between the incident flux h and the scattered flux I s : 

I s =(l-—)li. (45) 

In the following two subsections, we consider the cases where the polarization state of the incident wave is linear and circular, 
and compute the polarization states of the scattered mode. 



3.1 Linearly polarized incident wave 

Consider a linearly polarized electromagnetic wave (5i = 0, a, = 10 and Xi = 90°), which excites stimulated Raman scattering 
in the magnetospheric plasma. Using the plasma and magnetic field parameters discussed in the previous section, we computed 
the growth rate of stimulated Raman scattering. For fcAo <C 1, the instability becomes quite strong and the Landau damping 
of Langmuir wave is minimal. Hence, the stimulated Raman scattering is resonantly excited. Figure 4 shows the behaviour 
of linear (solid line) and circular (broken line) polarization of scattered mode with respect to rj. It shows at some values 
of rj, which are close to 0.1, the linear polarization of the incident wave can be converted almost completely into circular 
polarization of scattered wave. Charge density variations within the pulse window could enhance the conversion efficiency of 
linear to circular polarization to vary. 





8,=0, a,=10, 5 s =7i/2 
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Figure 4. The solid line and broken line curves indicate the variation of linear (L s ) and circular (V s ) polarization of the scattered mode 
with respect to rj. The normalizing parameter I s is the intensity of scattered mode. 



The variation of polarization angle of the scattered wave with respect to rj, at different values of 5, is shown in Fig. 5. For 
S = 90° and rj < 0.1, the scattered mode becomes orthogonally polarized with respect to the incident wave. Furthermore, if 
there is any variation in the plasma density or radiation intensity, the value of rj fluctuates and the scattered modes produced 
with rj < 0.1 become orthogonally polarized with respect to those produced with rj > 0.1. 
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Figure 5. Polarization angle \ s of the scattered mode is plotted as a function of rj, at different values of <5 (0°, 30°, 60° & 90°). 



3.2 Circularly polarized incident wave 

Suppose that the incident wave is circularly polarized (Si = tt/2 and q 4 = 1) then the scattered mode will be linearly polarized 
for t] < 0.2, while at the other values, both linear and circular polarizations with different proportion exist, as indicated by 
Fig. 6. For Si < and S s — Si — 8 < the sense of circular polarization of the scattered mode becomes opposite to that of 
the incident wave. Hence depending upon the plasma and radiation conditions, stimulated Raman scattering can change the 
polarization of the pulsar radio signals. 
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Figure 6. The variation of linear (L s ) and circular (V s ) polarization of the scattered mode, with respect to rj, arc indicated by the solid 
line and broken line curves. The normalizing parameter Is is the intensity of the scattered mode. 
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4 DISCUSSION 



The variable nature of circular polarization is evident from the polarization distributions by Manchester, Taylor and Huguenin 
(1975), Backer and Rankin (1980) and Stinebring et al. (1984a,b). Very high degrees of circular polarization are occasionally 
observed in individual pulses, even up to 100 per cent (Cognard et al. 1996), while the integrated or average pulse profiles 
generally indicate much smaller degree of circular polarization (e.g., Lyne, Smith & Graham 1971), which show that the sign 
of circular polarization is variable at any given pulse phase. Radhakrishnan and Rankin (1990) have identified two extreme 
types of circular polarization in the observations: (a) an antisymmetric type wherein the circular polarization changes sense 
in midpulse, and (b) a symmetric type wherein it is predominantly of one sense. The correlation of sense of antisymmetric 
type of circular polarization with the polarization angle swing indicates the geometric property of emission process, and is 
highly suggestive of curvature radiation. 

The diverse nature of circular polarization may be the consequence of pulsar emission mechanism and the subsequent 
propagation effects in the pulsar magnetosphere (e.g., Melrose 1995; von Hoensbroech, Lesch and Kunzl 1998). It seems 
rather difficult to explain the various circular polarization behaviours using the widely accepted magnetic pole models (e.g., 
Radhakrishnan & Cooke 1969; Komesaroff 1970; Sturrock 1971; Ruderman & Sutherland 1975). 

The symmetric type of circular polarization observed in some pulsars (e.g., PSR B1914+13, 0628-28) as shown in Figs. 7 
and 8, may be due to the propagation effects. However, it seems difficult for propagations effects to explain, how the sign of 
the circular polarization can change precisely at the center of the pulse in the case of antisymmetric type, as seen in many 
pulsars e.g., PSR B1859+03 and B1933+16 (Rankin, Stinebring & Weisberg 1989). 
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Figure 7. PSR B0628-28, an example of a pulsar with 'symmetric' circular polarization and high linear polarization (From Lyne & 
Manchester 1988). 

If we are to understand the radio emission mechanism, we must understand the physical state of the radio-loud plasma 
in the polar cap. It is this plasma that is the site of instabilities which are thought to produce coherent radio emission. The 
role of different propagation effects on the pulsar polarization has been discussed by Cheng and Ruderman (1979), Beskin, 
Gurevich and Istomin (1988) and Istomin (1992). The mechanisms proposed by these authors predict a frequency dependence 
for circular polarization, with weaker polarization at higher frequencies. This is seen in some pulsars (e.g., PSR B0835-41, 
1749-28, 1240-64) but it is not generally the case (Han et al. 1998). Istomin suggested that the linearly polarized incident 
wave becomes circularly polarized as a result of generalized Faraday rotation, however, it is observationally known that no 
generalized Faraday rotation is evident in pulsar magnetospheres (Cordes 1983; Lyne & Smith 1990). 
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Figure 8. Polarization of PSR B1914+13, a pulsar with strong circular polarization over the whole observed profile (From Rankin, 
Stinebring & Wcisbcrg 1989). 

We have presented a model to explain the polarization changes due to the propagation of radio waves through the 
magnetospheric plasma. The features like a large change in polarization angle, sense reversal of circular polarization and 
extremely rapid temporal changes in intensity would help us to explain many observations, for which, the existing mechanisms 
proved to be inadequate. Because of the very strong dependence of polarization angle on plasma parameters via the growth 
rate, in an inhomogeneous plasma medium the depolarization is a natural outcome. We believe that the plasma process such 
as the stimulated Raman scattering may be a potential mechanism for the polarization variability in pulsars. The circular 
polarization of a number of pulsars varies with frequency. The two clear examples are PSR B1240-64 and 2048-78, from which 
the opposite senses or transitions of circular polarization have been observed at lower and higher frequencies (Han et al. 1998). 



5 CONCLUSION 

We considered the stimulated Raman scattering of the transverse electromagnetic waves in the electron-positron plasma of 
pulsar magnetosphere. The value of radiation-plasma coupling parameter r\ is determined by the polarization of the incident 
wave, and its value can be determined only by the nonlinear analysis. In the process of three-wave interaction, the phase 
matching condition (equation 13) between the initial phases (5i, 8 S , S) and the value of rj determine the polarization state of 
the scattered wave. 

Many short time scale polarization variabilities, in individual pulses, such as variations in amount of linear and circular 
polarization, sense reversal of circular polarization and polarization angle swings can be accounted for by considering the 
stimulated Raman scattering. The time scales over which the changes takes place, in individual pulses, are of the order of 
e-folding time, which is the inverse of the growth rate of stimulated Raman scattering. 

It seems rather unlikely that the diverse behaviour of circular polarization can be accounted for by a single mechanism. 
Both intrinsic emission and propagation effects seem possible. The strong symmetric type of circular polarization observed in 
some pulsars is most probably generated by propagation effects, such as the stimulated Raman scattering. Further simultaneous 
observations over a wide frequency range would be valuable in sorting out the importance of propagation effects. 
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